Qubit-induced phonon blockade as a signature of quantum behavior in nanomechanical resonators 
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The observation of quantized nanomechanical oscillations by detecting femtometer-scale displacements is a 
significant challenge for experimentalists. We propose that phonon blockade can serve as a signature of quantum 
behavior in nanomechanical resonators. In analogy to photon blockade and Coulomb blockade for electrons, 
the main idea for phonon blockade is that the second phonon cannot be excited when there is one phonon in the 
nonlinear oscillator. To realize phonon blockade, a superconducting quantum two-level system is coupled to the 
nanomechanical resonator and is used to induce the phonon self-interaction. Using Monte Carlo simulations, the 
dynamics of the induced nonlinear oscillator is studied via the Cahill-Glauber s-parametrized quasiprobability 
distributions. We show how the oscillation of the resonator can occur in the quantum regime and demonstrate 
how the phonon blockade can be observed with currently accessible experimental parameters. 

PACS numbers: 85.85.+J, 03.65.Yz, 85.25.Cp, 42.50.Dv 



I. INTRODUCTION 

Many efforts (e.g., see Refs. 01-0] and reviews &0|) have 
been made to explore quantum effects in nanomechanical 
resonators (NAMRs) and optomechanical systems (e.g., in 
Refs. |8Hl3|] and the review AMI '). Reaching the quantum 
limit of NAMRs would have important applications in, e.g., 
small mass or weak-force detections fl5l - [l7ll . quantum mea- 
surements [18], and quantum-information processing. Only 
recently the quantum limit in NAMRs has been reached ex- 
perimentally 111911 . 

Quantum or classical behavior of a NAMR oscillation de- 
pends on its environment, which induces the decoherence and 
dissipation of the NAMR states. In principle, if the NAMR 
is cooled to very low temperatures (in the mK-range) and has 
sufficiently high oscillation frequencies (in the GHz-range), 
then its oscillation can approach the quantum limit. In other 
words, if the energies of the NAMR quanta, which are re- 
ferred to as phonons [20], are larger than (or at least compara- 
ble to) the thermal energy, then the mechanical oscillation can 
be regarded quantum. When the NAMR can beat the thermal 
energy and approach the quantum regime, measurements on 
quantum oscillation of the NAMRs are still very challenging. 
One encounters: (i) fundamental problems as measurements 
are usually performed by the position detection, the quan- 
tum uncertainty due to the zero-point fluctuation will limit 
the measurement accuracy; (ii) practical problems as, for a 
beam oscillating with frequency in the gigahertz range, the 
typical displacement for this oscillation is on the order of a 
femtometer. Detecting so tiny displacement is a difficult task 
for current experimental techniques. 

Various signatures and applications of quantum behavior 
(or nonclassicality) in nanomechanical resonators have been 
studied. Examples include: generation of quantum entangle- 



ment 1 2 1142311. generation of squeezed states Il24l - l27ll . Fock 
states ll28l l29ll . Schrodinger cat states l30ll . and other non- 
classical states 1I3TI l32tl . prediction of classical-like J33ll and 
quantum ll3~4ll Rabi oscillations, transport measurements 1 35 1 
quantum nondemolition measurements lfl8l l28l I29I l36i 
quantum tunneling 113 811 . proposal of quantum metrology 
and of quantum decoherence engineering ll25ll . 

The problem of how to perform quantum measurements 
on a system containing a NAMR plays a fundamental role in 
reaching the quantum limit of the NAMR and testing its non- 
classical behavior. Quantum measurements are usually done 
by coupling an external probe (detector) to the NAMR (see, 
e.g., IS 53-611] and references therein). 

Our approach for detecting quantum oscillations of 
NAMRs is based on: (i) recent theoretical proposals (e.g., 
Refs. (45J) to perform quantum measurement on NAMR with- 
out using an external probe and (ii) experimental demonstra- 
tions (e.g., Refs. ||46[ l47ll ) on the couplings between super- 
conducting quantum devices and the NAMRs. Instead of di- 
rectly detecting a tiny displacement, we propose to indirectly 
observe quantum oscillations of the NAMR via phonon block- 
ade, which is a purely quantum phenomenon. 

We assume that the phonon decay rate is much smaller than 
the phonon self-interaction strength. In such a case, we show 
that when the oscillations of the NAMR are in the nonclassical 
regime, the phonon excitation can be blockaded. In analogy 
to the photon (e.g., see Refs. I48ll49ll ) and Coulomb (e.g., see 
Ref. IdOiI ) blockades, the main idea for the phonon blockade 
is that the second phonon cannot be excited when there is one 
phonon in the nonlinear oscillator. Therefore, by analyzing 
correlation spectra for the electromotive force generated be- 
tween two ends of the NAMR, the phonon blockade can be 
distinguished from excitations of two or more phonons. 

An important ingredient for the realization of the phonon 
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FIG. 1 : (Color online) Schematic diagram for the coupling between 
a nanomechanical resonator (light blue bar on the right side) and su- 
perconducting charge qubit (left side). Two black dashed curves on 
the right show that the resonator is oscillating. The static magnetic 
field B (presented by the black upward-pointing arrow) and the al- 
ternating current I(t) (shown by the red loop on the right) are used 
for the motion detection of the NAMR. 



blockade is strong phonon self-interaction. To obtain such 
nonlinear phonon-phonon interaction, the NAMR is assumed 
to be coupled to a superconducting two-level system, which 
can be either charge, flux, or phase qubit circuits [51-54]. By 
choosing appropriate parameters of two-level systems, a non- 
linear phonon interaction can be induced. The interactions 
between each of these qubits [51-54] with NAMRs are very 
similar, e.g., the coupling constants are of the same order and 
the frequencies of these qubits are in the same GHz range. 
Therefore, in this paper, we only use charge qubits as an exam- 
ple to demonstrate our approach. However, this approach can 
also be applied to demonstrate the oscillation of the NAMR 
in the quantum regime, when the NAMR is coupled to other 
superconducting qubits, e.g., phase or flux qubits. 

The paper is organized as follows. In Sec. II, we describe 
the couplings between the superconducting qubits and the 
NAMR, and then study how the qubit induces the phonon- 
phonon interaction. In Sec. Ill, we discuss how to charac- 
terize the quantum oscillation by using the Cahill-Glauber s- 
parametrized quasiprobability distributions for s > 0, in con- 
trast to the Wigner function (for s — 0). In Sec. IV, the ba- 
sic principle of the phonon blockade is demonstrated, and we 
show that the phonon blockade can occur for the different pa- 
rameters. In Sec. V, we study the measurement of the phonon 
blockade by using the correlation spectrum of the electromo- 
tive force between two ends of the NAMR. Finally, we sum- 
mary the main results of the paper in Sec. VI. 



II. QUBIT-INDUCED PHONON-PHONON INTERACTION 

Let us now focus on the coupling between a NAMR (with 
mass TO and length L) and a superconducting charge qubit 
(with Josephson energy Ej and junction capacitance Cj). As 
schematically shown in Fig. Q] a direct-current (d.c.) voltage 
V g , and an a.c. voltage V g (t) — Vo cos(wit) are applied to the 
charge qubit (or Cooper pair box) through the gate capacitor 
C g . The NAMR is coupled to the charge qubit by applying a 



static voltage V x through the capacitor C(x) which depends 
on the displacement x of the NAMR around its equilibrium 
position. A weak detecting current I(t) — Iq cos(o>2i) is ap- 
plied to the NAMR, with its long axis perpendicular to the 
static magnetic field B. 

In the rotating wave approximation and neglecting two- 
phonon terms, the Hamiltonian H = Hq + of the interac- 
tion system between the charge qubit and the NAMR can be 
described by JH]: 

— — frcj a z + hcja f a + hg(aa + +a^cr_) 

+tin (o-+e- iUlt + a-e™ 1 *) , (1) 
7j( d ) = he (aV^ 2 * + ae iU2t ) . (2) 

Here, the frequency shift of the NAMR, due to its coupling to 
the charge qubit, has been neglected because it just renormal- 
izes the NAMR frequency and will not affect the calculations 
below. This frequency shift is determined lf3~lTl by the qubit- 
NAMR distance I, the charging energy E c = e 2 /2(Cj + C g + 
C), the mass m and the frequency cj of the NAMR. It should 
be noted that below we consider the large detuning between 
the qubit and the NAMR, i.e., (cjo — cj) is several times larger 
(but not much larger) than the coupling constant g; thus, the 
rotating wave approximation can be applied. The effect of 
the counter-rotating terms on the results can also calculated 
in the large detuning case 15511 . However, here we have ne- 
glected this effect because it only produces a small frequency 
shift and two-photon processes. The charge qubit, described 
by the spin operator a z = |e)(e| — |5)(<?|, is assumed to be 
near the optimal point, i.e., (C g V g + CV x )/2e ps 0.5 with 
C = C(x = 0), and thus cjq ~ Ej/K. The qubit ground and 
excited states are denoted by \g) and |e), respectively. The 
operator a (aJ) denotes the annihilation (creation) operator of 
the NAMR with frequency to, which can be written as 

Imcj ( i \ 

with the momentum operator p of the NAMR. The third 
term of Eq. (Q]l presents the NAMR-qubit interaction with the 
strength 

AE C N X X 



determined by the charging energy E c , effective Copper pair 
number N x = CV x /2e, the distance d between the NAMR 
and the superconducting qubit, and the NAMR amplitude 
Xa = y/h/2mcj of zero-point motion. Also, ft is the Rabi 
frequency of the qubit driven by the classical field with fre- 
quency cji . The parameter 

e = -BI LX (6) 

in Eq. $2^ describes the interaction strength between the 
NAMR and an external weak probe a.c. current with fre- 
quency cj2- Hereafter, we assume that the resonant driving 
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condition for the qubit is satisfied, i.e., u>y = loq. For the cou- 
pling between a phase [ 2 1 ] (or flux qubit 1 46 ] ) and the NAMR, 
they also have the same form as that given in Eqs. ([T]i and iff} 
except all parameters of the Hamiltonian should be specified 
to the concrete systems. Thus, our discussions below can also 
be applied to those systems. 

The frequency of the NAMR is usually much lower than 
that of the qubit. If the Rabi frequency satisfies the condi- 
tion ft ^> (g 2 /A) with the detuning A = luq — cj between 
the frequencies of the NAMR and the qubit, then in the ro- 
tating reference frame with V = exp(—iuioa z /2), Eq. (Q]) is 
equivalent to an effective Hamiltonian 



H. 



(0) 



ctr 



9 t i t \i 
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(7) 



with the effective phonon self-interaction constant (nonlinear- 
ity constant) 



QA 2 



(8) 



Here, p z 



-)(- 



-)(— | with the dressed qubit states 



|±) = (\g) ± |e) ) /a/2- Therefore, if the dressed charge qubit, 
which was theoretically proposed [56] and has been experi- 
mentally realized lf57l 15811 . is always in its ground state |— ), 
the effective Hamiltonian for the driven NAMR is 

H c ff = H ( w — ) a^a + hK(a^a) 



7ie(aV ia;2i -r-ae ia,2i ) 



(9) 

The nonlinear Hamiltonian of the driven NAMR in Eq. (0 
can also be directly obtained when the driving field is strong; 
however, here we only consider a weak probe current. Thus, 
the coupling of the NAMR to a controllable superconducting 
two-level is necessary for inducing phonon-phonon interac- 
tions. 



HI. QUANTUM BEHAVIOR DESCRIBED BY 
QUASIPROBABILITY DISTRIBUTIONS 

Decoherence imposes strict conditions to observe quantum 
behavior in a NAMR. To demonstrate effects of the environ- 
mental on the NAMR, let us now assume that the NAMR is 
coupled to the thermal reservoir. Under the Markov approxi- 
mation, the evolution of the reduced density operator p for the 
NAMR can be described by the master equation [59]: 



d_ 

dt 



p 



T - t 
:[H cS , p] + -n(2a T pa 



aa 1 p - 



paa 



l)(2apa f — aJap — pa' a) 



(10) 



In Eq. ([Tol l, 7 is the damping rate and n = {exp[fiw / '{ksT)\ — 
is the mean number of thermal phonons, where ks is 
the Boltzmann constant, and T is the reservoir temperature at 
thermal equilibrium. Eq. (TlOb can be solved , e.g., b y applying 
the Monte Carlo wave function simulation 0591 - 16111 and intro- 
ducing the collapse operators 



d = ^j{n + l)a , C 2 



(a) 




X 




•fna' 



(11) 
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FIG. 2: (Color online) Quasidistribution functions for the NAMR 
steady state obtained by solving master equation J 10b for e = 
37, k = 3O7, and n = 0.01 with 7 as units: (a) Wigner 
function W^°\x,y), which is non-negative in the whole phase 
space, and (b) 1/2-parametrized quasi-probability distribution (QPD) 
VF' 1 / 2 ' (x, y), which is negative for a = x + iy close to zero indicat- 
ing nonclassicality of the NAMR state. The figures show the bottom 
of the functions. 



We now study the steady-state solution, which is indepen- 
dent of the initial states. For the system without a drive, the 
time evolution ends in a state without phonons (vacuum state) 
at zero temperature. While for a driven system, the asymp- 
totic state is neither the vacuum nor a pure state even at zero 
temperature, and can have intriguing noise properties. 

A state is considered to be nonclassical if its Glauber- 
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Sudarshan P function cannot be interpreted as a probability 
density, i.e., it is negative or more singular than Dirac's S 
function. Due to such singularities, it is usually hard to vi- 
sualize it. To characterize the nonclassical behavior of the 
NAMR states generated in our system, we consider the Cahill- 
Glauber s-parametrized quasiprobability distribution (QPD) 
functions 16211: 



= -Tc\pT^(a)] 



(12) 



where 



= \J expK* -a*0^ (s) (Od 2 C, (13) 



and 



with 



^IO=«*i>(*1^)2>(0. 



D{g) = exp _ £* a ) 



(14) 



(15) 



being the displacement operator. The QPD is defined for 
— 1 < s < 1, which in special cases reduces to the P function 
(for s = 1), Wigner function (for s = 0), and Husimi Q func- 
tion (for s — —1). QPDs contain the full information about 
states. 

Let us analyze the differences between the 1/2- 
parametrized QPD and Wigner functions under the resonant 
driving for the NAMR with bj 2 = lj - (g 2 /A). As an 
example, in Fig. [2] we plotted the steady-state Wigner 
function and 1/2 -parametrized QPD, which are the numerical 
solutions of the master equation for a set of parameters: 
n = 0.01, e = 37, and k — 3O7 in units of 7. Fig. |2|a) 
shows the non-negative Wigner function of the steady state 
of the NAMR for these parameters. It can also be shown 
analytically that the steady-state Wigner function for this 
system is always non-negative. However, the plot for the QPD 
function W (1/2) (a) in Fig. Eft), with the same parameters 
as for Fig.|2a), clearly shows negative values, corresponding 
to a nonclassical state of the NAMR. Below, we will discuss 
how to demonstrate this nonclassicality of the NAMR via the 
phonon blockade. 

The Wigner function for the NAMR steady state is non- 
negative in the whole phase space. This is in contrast to 
the Wigner function for various nonclassical states, including 
Fock states or finite superpositions of coherent states (often 
referred to as Schrodinger cat states) being negative in some 
areas of phase space. It should be noted that there are other 
well-known nonclassical states, including squeezed states, for 
which the Wigner function is non-negative as for the NAMR 
steady state. In general, the non-positivity of the Wigner 
function is a necessary but not a sufficient condition for the 
non-classicality. The complete characterization of the non- 
classicality (the "if and only if" condition) is based on the 
positivity of the P-function. Unfortunately, this function is 
usually too singular to be presented graphically. The larger 
parameter s the more nonclassical states are described by the 




FIG. 3: (Color online) Probabilities P n = (n\p{i)\n) of measuring 
n phonons as a function of rescaled time, et, assuming k = lOe 
and: (a) no dissipation (7 = 0) and (b) including dissipation with 
the same parameters as in Fig. [2] Po (red curves), Pi (blue), and 
P2 (green). F — Po + Pi (thick black) describes the fidelity of 
the phonon blockade. Additionally, the coherences X — Re(0|p|l) 
(magenta curves) and Y = lm(0|p|l) (cyan) show that the steady 
states partially preserve coherence. 



negative s-parametrized QPD. In our case, to demonstrate the 
nonclassically of the NAMR steady state, it was enough to 
calculate the s-parametrized QPD for s = 1/2 but not for 
s = 0. 



IV. PHONON BLOCKADE 

We now consider the case when the phonon self-interaction 
strength k is much larger than the phonon decay rate 7. 
When the oscillation of the NAMR is in the quantum regime, 
the phonon transmission can be blockaded in analogy to the 
single-photon blockade in a cavity [48, 49]. This is because 
the existence of the second phonon requires an additional en- 
ergy Hk. 

To demonstrate the phonon blockade, let us rewrite Eq. (O 
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H eS = huja)a + hKa 1 a(<Ja - 1) + 7ie(a t e" ia ' 2t + ae iui2t ) 

(16) 

with a renormalized frequency 



UJ = UJ + K — — . 



(17) 



In the rotating reference frame for V' — exp(— iwza'a t) with 
uj'2 = u, the Hamiltonian in Eq. (fToT ) becomes 



H e s = Hna) a(a' a — 1) + he(a ! + i 



(18) 



It is now easy to see that the two states |0) and |1) with zero 
eigenvalues are degenerate in the first term na^a^a — 1) of 
Eq. (1181 . This degeneracy plays a crucial role in the phonon 
blockade. Indeed, if we assume that the interaction strength e 
is much smaller than the nonlinearity constant k (i.e., e <C k), 
then the phonon eigenstates of the Hamiltonian in Eq. ( fT8l can 
become a superposition of only two states, |0) and |1), in the 
lowest-order approximation of the expansion in the strength e. 

We now study the solution of the Hamiltonian in Eq. ( TToT l 
under the assumption of a weak driving current, i.e., e <C K. 
Using standard perturbation theory, the state governed by the 
time-dependent periodic Hamiltonian in Eq. ( TToT l with the ini- 
tial condition \ip(t = 0)) = |0) can be obtained by introducing 
the auxiliary operator 



d 

H F = H ca - i — 
at 



(19) 



based on the Floquet theory (e.g., see Ref. 116311 ). The solution 
can be approximately given as 

|t/>(t)) =cos(et)|0) -£sin(ei)|l) +C(e 2 ). (20) 

The solution in Eq. (1201 shows that the number of phonons 
varies between and 1 if all terms proportional to e 2 are ne- 
glected. In this small e limit, the Floquet solution (|20| | explic- 
itly demonstrates the phonon blockade in analogy to the pho- 
ton blockade J49ll or the Coulomb blockade l50ll . i.e., there is 
only one-phonon excitation and the excitation with more than 
one phonon is negligibly small. The photon blockade is also 
referred to as the optical state truncation 11481 164I1 . 

The phonon-blockaded state is nonclassical as it is a super- 
position of a finite number (practically two) of Fock states. 
Only (some) superpositions of an infinite number of Fock 
states can be considered classical. 

The time-dependent probabilities P n = (n\p(t)\n) of mea- 
suring the n-phonon state with and without dissipation are nu- 
merically simulated using the Monte Carlo method. In the 
ideal non-dissipative case, as shown in Fig. 0a), the sum of 
the probabilities Po and P\ with phonon numbers and 1 is 
almost one, which means that phonon blockade occurs. For 
the dissipative case, Fig.[3jb) shows the time evolutions of the 
elements (m\p(t)\n) (with m, n — 0, 1) for the same param- 
eters as those in Fig. [2] The amplitudes of P an d Pi exhibit 
decaying oscillations; however, their sum is still near one and 
thus the sum of other probabilities P n with n > 1 is near zero. 
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FIG. 4: (Color online) Probabilities P n , fidelity F, and coherences 
X and Y for steady states as a function of (a) /3 = huj / '(ksT), 
assuming k — We, and of (b) ft/7, assuming n = 0.01, which 
corresponds to j3 f» 4.6. In both (a) and (b), we set 7=1. The other 
parameters are the same as in Fig. [3] 



Therefore, the phonon blockade occurs even in the long-time 
limit (e.g., steady state). The non-zero off-diagonal element 
(0|jo(f)|l) shown in Fig. Ob) in the steady-state means that 
the NAMR is in the nonclassical state, which is also consis- 
tent with the steady-state plot of the QPD in Fig.|2|b). Thus, 
we see that the non-negative Wigner function does not directly 
indicate that the state is nonclassical. 

To study how the environmental temperature T affects 
the phonon blockade, the probability distributions P n = 
(n\p(t)\n) (for n = 0, 1, 2, 3) are plotted via the rescaled 
inverted temperature Sa;/(fcsT) in Fig. |4f a). It clearly shows 
that the phonon blockade cannot be achieved when the ther- 
mal energy is much larger than that of the NAMR. The k- 
dependent matrix elements (m\p(t)\n) are plotted in Fig. SJb), 
which shows that the larger nonlinearity parameter n corre- 
sponds to a more effective phonon blockade. However, to ob- 
serve the phonon blockade, it is enough to make k larger than 
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a certain value. For instance, if the ratio between k/j is larger 
than 10, then the sum of the probabilities P and Pi is more 
than 0.95, and the phonon blockade should occur. 

Let us make a few comments to clarify the relation between 
the phonon blockade and nonclassicality in terms of the s- 
parametrized QPDs: (i) If the s-parametrized QPD, for some 
s € (—1, 1] and for a given state, is negative in some region 
of the phase space, then the state is nonclassical. (ii) Even 
if the phonon blockade is not observed (for a given choice of 
parameters e, k, 7, and n), the 1/2-parametrized QPD (or the 
QPD for any s > — 1) can still be nonpositive. (iii) Even if we 
choose the parameters such that the 1/2-parametrized QPD is 
positive, this does not imply that the state is classical, (iv) 
Even if the phonon blockade does not appear, the state can be 
nonclassical as described by the nonpositive P-function (the 
QPD for s = 1). 

A good blockade of phonons can be observed for nonclassi- 
cal states only. However, a poor blockade of phonons does not 
imply that the state is classical. Similarly to other quantum ef- 
fects like squeezing or antibunching: If a specific nonclassical 
effect is not exhibited by a given state, it does not imply that 
the state is classical. 

We can choose the parameters e, k, 7, and n in order 
to observe a change (transition) from a nonpositive 1/2- 
parametrized QPD to a positive function. However, this tran- 
sition is not important in the context of nonclassicality. For 
various s > —1, one could observe such transitions for dif- 
ferent parameters. Only the transition of the P-function cor- 
responds to a transition from quantum to classical regime. As 
already mentioned, a good criterion of nonclassicality should 
be based on the P-function, but it is usually too singular to be 
presented graphically. Thus, we have chosen the QPD for an- 
other value of s € (0, 1). A nonclassicality criterion based on 
the QPD for s = 1/2 is more sensitive than that based on the 
Wigner function (the QPD for s = 0), but still it is not sen- 
sitive enough in the general case, i.e., there are nonclassical 
fields described by the positive 1/2-parametrized QPD. 



V. PROPOSED MEASUREMENTS OF THE PHONON 
BLOCKADE 

Let us now discuss how to measure the phonon blockade 
via the magnetomotive technique, which is one of the basic 
methods to detect the motion of NAMRs l65ll . As shown in 
Fig. Q] the induced electromotive force V between two ends 
oftheNAMR is HHH 

V = BL^- = %BL\\ — (a* - a) , (21) 
m V 2m 

which can be experimentally measured as discussed in 
Ref. l65ll . Here, p is the momentum for the center of the 
NAMR mass. We analyze the power spectrum 

/CO 
{V{Q)V{T))e- lbJ,t dt (22) 
-00 



defined by the Fourier transform of the induced electromotive- 
force two-time correlation function 

(V(O)V(t)) = lim (V(t)V(t + r)) . (23) 
t— ¥00 

This power spectrum can be measured effectively. 

Power spectrum Sv(u') and the two-time correlation func- 
tion (V(0)V(t)) are plotted for zero temperature with differ- 
ent decay rates 7 in Fig. |3a), and for a given decay rate with 
different temperatures T (i.e., different thermal phonon num- 
ber n) in Fig. [2b). We find that low dissipation and low tem- 
peratures produce high spectral peaks, which enable an easier 
observation of the phonon blockade. Thus, the environment 
(or some background) will limit the power spectrum for ob- 
serving the phonon blockade. When k is negligible compared 
with the decay rate 7, all spectral peaks disappear and there is 
no phonon blockade. By other numerical calculations, we also 
find that a large or giant nonlinearity k corresponds to sharp 
peaks and, in this case, the phonon blockade is also easy to be 
observed. 

Assuming perfect phonon blockade, i.e., truncation to an 
exact qubit state, one can analyze the whole evolution of our 
system confined in two-dimensional Hilbert space. To some 
extend this approximation can be applied in our model if the 
conditions e k and (n) w are satisfied. Then, we find that 
the corresponding power spectrum should have at most three 
peaks at frequencies 

u' = 0, uo' 12 = ±iy / (8e) 2 -7 2 (l + 2n) 2 « ±2e. (24) 

It is seen that these frequencies are independent of k. A peak 
at uj' q = does not appear for real e, which is the case ana- 
lyzed in the paper. Examples of such power spectra for u/ > 
are shown in Fig. 5(a) for n = and in Fig. 5(b) for n = 0.01 
(blue curve). In contrast, new peaks appear in the spectrum 
in the case of not perfect phonon blockade. This can be un- 
derstood by analyzing a Hilbert space of dimension d > 2. 
For example, by analyzing the system evolution confined in 
a three-dimensional Hilbert space, we find that the spectrum 
can have at most seven peaks centered at 

J k. 0, ±2e (1 - 5) , ± [2k (1 + 65) ± e (1 - 5)] , (25) 

where S = 6 2 /(8k 2 ), which depend on n, contrary to 
Eq. (l24l . Frequencies in Eq. ( l25l > can be approximated as 
u)' « 0,±2e,±(2K ± e). Thus, for lu' > 0, the first peak 
occurs at 2e, which corresponds approximately to uu[ given in 
Eq. d24l >. The second characteristic double peak is at 2k± e, as 
seen in Fig. 5(b) for n — 0.5 (red) and n = 1 (black curves). 
Eq. (|25l l explains only the occurrence of the first three peaks 
for uj' > Q in Fig. 5(b). To explain the appearance of the other 
two peaks at lu 1 k, 4k and 6k, one should analyze the evo- 
lution of our system confined in (at least) four-dimensional 
Hilbert space. Thus, these extra peaks are a signature of a 
non-perfect single phonon blockade. 

The spectra are not symmetric in frequency around zero, 
SV(u/) 7^ Sv{— w'). Nevertheless, we depicted only the 
positive-frequency half of the spectra in Fig. 5 to better com- 
pare the peaks for different values of n. We note that a double 
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peak is observed at negative frequencies u/ — (2k ± e) even 
for the cases shown in Fig. 5(a). This means that the contribu- 
tion of terms 0(e 2 ) in Eq. ( f20l > is not negligible for the param- 
eters chosen in Fig. 5(a), and the spectrum for oj' < does 
not correspond to a (mathematically) perfect single-phonon 
blockade. 

In Fig. 5(b), the power spectra are plotted as a function of 
u/ j k. There, it is seen that the position of the first positive 
peak depends on the ratio e/n in agreement with Eq. ( f25l ). 
The center of this peak is closer to zero for smaller ratio e/n. 
However, the position of the center of the double peak (split 
peak) is approximately independent of e and k (assuming that 
k e, so 5 « 0), which follows from Eq. (1251 , Moreover, the 
splitting vanishes with increasing 7. The smaller e the smaller 
is 7 for which the splitting vanishes. 

In conclusion, the observation of extra peaks at frequen- 
cies different from those in Eq. (l24t . show deterioration of the 
single-phonon blockade. The higher are such peaks the worse 
is the phonon blockade. 

Note that the double peak at uj' w 2k ± e was found as- 
suming the output state to be in a qutrit (three-dimensional) 
state. This double peak can, in general, be predicted for a 
qudit state, i.e., e?-dimensional state for 2 < d <C 00. This 
corresponds to a phonon-truncation up to state \d — 1) and can 
be interpreted as a generalized multi-phonon blockade. Any 
qudit states are nonclassical since arbitrary finite superposi- 
tions of number states are nonclassical. However, with in- 
creasing dimension d of qudit states it becomes more difficult 
to distinguish them from classical infinite-dimensional states 
generated in our system. For this reason, here we analyze the 
standard single-phonon blockade only. 

We now discuss the experimental feasibility of our pro- 
posal. With current experiments on coupling a superconduct- 
ing phase JH (or charge JHHHHl] or flux Jil) qubit and 
the NAMR, the coupling constants are on the order of hun- 
dreds of MHz (e.g., 200 MHz), the environmental tempera- 
ture can reach several tens of mili-Kelvin (e.g., 20 mK), the 
frequency of the NAMR can be in the range of GHz (e.g., 1 
GHz). If the qubit frequency luq and the Rabi frequency fl 
are taken as, e.g., ujo = 2 GHz and il = 200 MHz, then 
the nonlinear parameter is k = 8 MHz. The observation of 
the phonon blockade should be possible for a quality factor Q 
larger than 10 3 , which is in the NAMR quality factor range 
10 3 ~ 10 6 of current experiments. By engineering k as in 
Refs. 1 66 68 1 . k can be much larger than 8 MHz, then the 
phonon blockade should be easier to be observed in our pro- 
posed system. 



VI. CONCLUSIONS 

We have studied the quantum mechanics of the NAMR by 
coupling it to a superconducting two-level system. To demon- 
strate our approach, a classical driving microwave is applied 
to the qubit so that a dressed qubit is formed. If the Rabi fre- 
quency of the driving field is strong enough, then the nonlinear 
phonon interaction can be induced when the dressed qubit is in 
its ground state. We mention that dressed charge qubits have 




FIG. 5: (Color online) Power spectra Sv(u') for k = 30, e = 3 and: 
(a) n = with 7 = 0.5 (blue), 1 (red), 1.5 (black), and (b) 7 = 0.5 
with n = 0.01 (blue), 0.5 (red), 1 (black curves). Parameters k, e 
and 7 are in units of g on the order of MHz. 



been experimentally realized ll57il58ll . The dressed phase (see, 
e.g., Refs. JUHzUl) and flux (see, e.g., Ref. SB) qubits 
should also be experimentally realizable. 

The states of the nonlinear NAMR can be completely char- 
acterized by the Cahill-Glauber s-parametrized quasiprobabil- 
ity distribution (QPD). A state is considered to be nonclassical 
if it is described by a P-function (QPD for s = 1) that can- 
not be interpreted as a probability density. As a drawback, the 
P-function is usually too singular to be presented graphically. 
Thus, other QPDs are often analyzed: If, for a given state, 
a QPD with s > —1 is negative in some regions of phase- 
space, then the state is nonclassical. We have shown that the 
Wigner function (QPD for s = 0) is always non-negative for 
nonclassical steady states generated in our dissipative system. 
Thus, we have calculated the 1/2-parametrized QPD, being 
negative in some regions of phase-space, which clearly in- 
dicates the nonclassical character of the steady states gener- 
ated in our NAMR system. Nevertheless, from an experimen- 
tal point of view, the quantum-state tomography of the 1/2- 
parametrized QPD is very challenging. Thus, we have pro- 
posed another experimentally-feasible test of nonclassicality: 
the phonon blockade. 

We considered the case when the phonon self-interaction 
strength k significantly exceeds the phonon decay rate 7. We 
showed that when the NAMR oscillations are in the quantum 
regime, the phonon transmission can be blockaded in analogy 
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to the single-photon blockade in a cavity [48, 49] or Coulomb 
blockade for electrons 115011 - When the phonon blockade hap- 
pens we also showed that a NAMR is in a nonclassical state 
even if its Wigner function is non-negative. Therefore, the 
nonclassicality of the NAMR can be demonstrated by the 
phonon blockade, instead of trying to detect the tiny displace- 
ments when the NAMR approaches the quantum limit. We 
further demonstrated that the phonon blockade can be exper- 
imentally observed by measuring the correlation spectrum of 
the electromotive force. All parameters in our approach are 
within current experimental regimes and, therefore, the quan- 
tum signature of the NAMR might be demonstrated in the near 
future by using this proposed approach. 

We have shown that the phonon blockade can be demon- 
strated by a qubit-induced nonlinear NAMR. However, the 
temperature of the environment, the decay rate of the NAMR, 
the driving current, and the nonlinear coupling constant k 
limit the measured power spectrum. To more efficiently ob- 
serve the phonon blockade, the following conditions should 
also be satisfied: (i) the temperature should be low enough so 
that thermal excitations should be negligibly small or the ther- 
mal energy is smaller than that of the oscillating energy of the 
NAMR; (ii) the quality factor of the NAMR should be high; 
(iii) the driving current through the NAMR should be very 
weak, so that the heating effect induced by the driving current 
can be neglected; (iv) the giant nonlinear constant k of the 
NAMR might be more useful for the phonon blockade, and 



this might be obtained using the approaches explored, e.g., in 
Refs. [ 66-68] . In our proposal, the larger coupling constant 
g between the qubit and the NAMR corresponds to a larger 
k, and the phonon blockade should be more easily observable 
for larger k. Also the frequency of the NAMR should be large 
enough, so that the qubit and the NAMR are in the large de- 
tuning regime, but the detuning should not be extremely large. 
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